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We present preliminary results using overlap fermions for the charmonium spectrum, in particular for hyperfine 
splitting. Simulations are performed on 16'^ x 72 lattices, with Wilson gauge action at (3 — 6.3345. Depending on 
how the scale is set, we obtain 104(5) MeV (using IP — IS) or 88(4) MeV (using ro=0.5 fm) for the hyperfine 
splitting. 

1. Introduction 

Overlap fermions have the following desirable 
features: 

• Exact chiral symmetry on the lattice. 

• No additive quark mass renormalization. 

• No flavor symmetry breaking. 

• No 0{a) error. 



• The O{m'^o?) and ©(AqcD^^a^) errors are 
also small, from dispersion relation and 
renormalization constants. 
The first two features are especially significant 
for light quarks. Many exciting results at low 
quark masses have been reported using overlap 
fermions. Here we want to make the point that 
using overlap fermions can also alleviate some 
problems with simulating heavy quarks. The last 
feature, demonstrated in |l], is an unexpected 
bonus in this regard. The key observation is that 
the discretization errors are only about 5% all the 
way upto ma 0.5. 

We make our case using hyperfine splitting in 
the charmonium system. It is known that with 
staggered quarks, there is an ambiguity about 
Nambu-Goldstone (NG) and non-NG mode for 
the ?7c, resulting in widely different estimates of 
hyperfine splitting - 51(6) MeV (non-NG) and 
404(4) MeV (NG) 0. NRQCD converges only 
slowly for charm Including 0{v^) terms 

changed the result from 96(2) MeV to 55(5) MeV. 



Wilson fermions have 0{a) errors. Hyperfine 
splitting is very sensitive to the coefficient of the 
correction term, csw- There are many studies 
0] using Wilson type valence quarks, including 
some with non-perturbative csw- The quenched 
clover estimate of hyperfine splitting has stabi- 
lized around 80 MeV. Result from a 2-f 1 dynam- 
ical simulation using tree-level csw still falls short 
of the experimental value by about 20% 0. 

Although costly, overlap fermions offer the best 
solution — they do not have flavor symmetry 
breaking, they guarantee that the 0{a) error is 
absent without any tuning and they are rela- 
tivistic. Furthermore, O(a^) errors appear to be 
small 

2. Simulation Details 

Our simulations are performed on 16'^ x 72 
isotropic lattices. We present results on 100 con- 
figurations. The Wilson gauge action is used at /3 
= 6.3345. For quenched configurations, the scales 
set from different physical quantities can differ 
considerably. Hyperfine splitting is very sensitive 
to the scale, so we report results for scale set us- 
ing two different quantities: and the (IP— 15) 
mass splitting in charmonium. From (IP— IS) we 
obtain a = 0.0501 fm. Using tq = 0.5 fm the scale 
is a = 0.0560 fm Pj. We use a multimass inverter 
to obtain propagators for 26 masses ranging from 
0.020-0.85 in lattice units. The bare masses cor- 
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Figure 1. Effect of inner loop precision on pseu- 
doscalar propagators for heavy quarks. We study 
output of one spin and one color for a single con- 
figuration for this illustration. Curves are slightly 
shifted for clarity. 




0.038 
0.036 
0.034 
0.032 
0.03 
0.028 
0.026 
0.024 



Effectile mass for rktio of JW td ri^, ma = 0.i5 



********* T I- 



10 



15 



20 



25 



30 



35 



Figure 2. Effective hyperfine splitting (mj/^ — 
rriri^) from the ratio of the vector to pseudoscalar 
correlators for two values of ma which bracket the 
charm mass. 

respond to 70 MeV - 3 GeV. The charm mass in 
lattice units is 0.365(1) for the (IF — IS) scale, 
and 0.431(1) for rp scale. Near these values of 
ma our discretization errors are estimated to be 
under 5%. 

For heavy quarks, mesonic two-point functions 
fall through many orders of magnitude, and be- 
come very small at the center of the lattice. An 
imprecise quark propagator results in drastically 
different values for these small numbers. Our ob- 
servation is that a good inner loop precision is 
crucial, as shown in Fig. To confirm that our 
outer loop precision is adequate, we increased it 
by an order of magnitude. This lead to less than 
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Figure 3. We fit the meson masses linearly in 
quark mass. Fits are shown for rjcJ/ip and he 
masses. 

a percent deviation in the (single configuration, 
single spin-color) pseudoscalar correlator. For our 
production runs, we choose an inner loop preci- 
sion of 10^* and outer loop precision of 10~^, 
with eigenvalue projection precision of 10"^''. At 
charm mass, the residual is about 2 x 10^^. 

3. Analysis 

Fig. 121 shows the effective hyperfine splitting 
from the ratio of vector to pseudoscalar corre- 
lator. We use two ways to set the scale - from 
the (IP — IS*) splitting in the charmonium sys- 
tem and from rg (using 0.5 fm). The experi- 
mental mjf^ is used to set m^ (in lattice units). 
This is straight-forward when the scale from vq 
is used. The singlet P mass mu^ is used for P, 
and (3mj/^ -I- mrj^)/A for S mass. The spin av- 
eraged (IP — IS) sphtting is used because it is 
expected to be insensitive to lattice artifacts. In 
this case, the determination of a and mca is en- 
tangled. The procedure we follow to disentangle 
these is as follows. As shown in Fig. O all hadron 
masses in lattice units arc fitted to a straight line, 
m^a — Afi-ma + Bh- Lattice spacing a and bare 
charm quark mass mca are two unknowns; mj/^ 
and to(1P — IS") in physical units are the two 
inputs. We solve for a and to obtain the 
values quoted earlier. If we were to use this a 
to get the Sommer scale, we would obtain rg = 
0.45 fm. Fig. 01 shows the charmonium spectrum 
in physical units. Interpolation for the hyperfine 
splitting is shown in Fig. [SI The fit form used is 
(mj — m,j^)a = Aj ^ma + B /ma. The intcrpo- 
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lation for rric for the two scales is shown. Note, 
hyperfine spUtting is higher for (IP —IS) scale in 
lattice units and also a is smaller, so the result in 
MeV is considerably higher for that case. Finally 
we summarize the results in the table below. All 
masses are in GeV. 
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Table 1 

Charmonium spectrum (GeV) 
4. Summary 

We have presented the first study of the char- 
monium spectrum using overlap fermions. We 
get a better agreement with the experimental 
spectrum using IP — IS scale rather than the 
ro scale. Our value for hyperfine splitting is 
104(5) MeV and 88(4) MeV using IP - IS and 
ro scale respectively. This is considerably higher 
than the quenched clover results. We make two 
approximations: quenching and excluding OZI 
suppressed diagrams. At the charm mass the sin- 
glet contribution appears to be small 0, though 
lattice calculations with smaller statistical and 
systematic errors are needed to settle this issue. 
If the unquenched clover simulations are any indi- 
cation, then it is possible that dynamical overlap 
with only connected insertion will reproduce the 
experimental result. 
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Figure 4. Charmonium spectrum in physical 
units. Results from both ro and IP — IS scales 
are shown. Note, for the latter scale, a linear 
combination of he and rjc masses, along with the 
J/ip mass, is used for input. 
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Figure 5. Hyperfine splitting as a function of 
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